Chapter 23 



Why does it work? 



Bloch; "Space is the field of linear operators." 
Heisenberg: "Nonsense, space is blue and birds fly 
through it." 

— Felix Bloch, Heisenberg and the early days of 
quantum mechanics 

(R. Artuso, H.H. Rugh and R Cvitanovic) 

As WE SHALL SEE, the tracB formulas and spectral determinants work well, 
sometimes very well. The question is: Why? And it still is. The heuris- 
tic manipulations of chapter 18 were naive and reckless, as we are facing 
infinite-dimensional vector spaces and singular integral kernels. 

We now outline the key ingredients of proofs that put the trace and determi- 
nant formulas on solid footing. This requires taking a closer look at the evolution 
operators from a mathematical point of view, since up to now we have talked 
about eigenvalues without any reference to what kind of a function space the cor- 
responding eigenfunctions belong to. We shall restrict our considerations to the 
spectral properties of the Penon-Frobenius operator for maps, as proofs for more 
general evolution operators follow along the same lines. What we refer to as a "the 
set of eigenvalues" acquires meaning only within a precisely specified functional 
setting: this sets the stage for a discussion of the analyticity properties of spectral 
determinants. In example 23.1 we compute explicitly the eigenspectrum for the 
three analytically tractable piecewise linear examples. In sect. 23.3 we review the 
basic facts of the classical Fredholm theory of integral equations. The program 
is sketched in sect. 23.4, motivated by an explicit study of eigenspectrum of the 
Bernoulli shift map, and in sect. 23.5 generalized to piecewise real-analytic hy- 
perbolic maps acting on appropriate densities. We show on a very simple example 
that the spectrum is quite sensitive to the regularity properties of the functions 
considered. 

For expanding and hyperbolic finite-subshift maps analyticity leads to a very 
strong result; not only do the determinants have better analyticity properties than 



452 



CHAPTER 23. WHY DOES IT WORK? 453 

the trace formulas, but the spectral determinants are singled out as entire functions 

in the complex s plane. remark 23.1 

The goal of this chapter is not to provide an exhaustive review of the rigorous the- 
ory of the Perron-Frobenius operators and their spectral determinants, but rather 
to give you a feeling for how our heuristic considerations can be put on a firm 
basis. The mathematics underpinning the theory is both hard and profound. 

If you ai^e primarily interested in applications of the periodic orbit theory, you 
should skip this chapter on the first reading. 



fast track: 
chapter 13, p. 266 



23.1 Linear maps: exact spectra 

We start gently; in example 23.1 we work out the exact eigenvalues and eigen- 
functions of the Perron-Frobenius operator for the simplest example of unstable, 
expanding dynamics, a linear 1 -dimensional map with one unstable fixed point. 
Ref. [23.6] shows that this can be carried over to (i-dimensions. Not only that, 
but in example 23.5 we compute the exact spectrum for the simplest example of a 
dynamical system with an infinity of unstable periodic orbits, the Bernoulli shift. 

Example 23.1 The simplest eigenspectrum - a single fixed point: In order to get 
some feeling for the determinants defined so formally in sect. 19.2, let us work out a 
trivial example: a repeller with only one expanding linear branch 

fix) ^ Ax, |A|>1, 

and only one fixed point x,, = 0. The action of the Perron-Frobenius operator (16.10) is 



£<P(y) = jdxSiy- Ax)(/>ix) = ^^O'/A) • (23.1) 



1 



From this one immediately gets that the monomials y are eigenfunctions: 

= /, Ic^ 0,1,2,... (23.2) 
IMA*:-" 



What are these eigenfunctions? Think of eigenfunctions of the Schrodinger 
equation: k labels the ^th eigenfunction x'^ in the same spirit in which the number 
of nodes labels the kth quantum-mechanical eigenfunction. A quantum-mechanical 
amplitude with more nodes has more variability, hence a higher kinetic energy. 
Analogously, for a Pen^on-Frobenius operator, a higher k eigenvalue l/IAIA*^ is 
getting exponentially smaller because densities that vary more rapidly decay more 
rapidly under the expanding action of the map. 
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Example 23.2 The trace formula for a single fixed point: The eigenvalues A * ' 
fall off exponentially with k, so the trace of £ is a convergent sum 



1 " 
lAI 



-k_ 1 _ 1 



|A|^ lAI(l-A-i) 1/(0)' 

in agreement with (18. 7). A similar result follows for powers of £, yielding the single- 
fixed point version of the trace formula for maps (18.10): 

y^^^y^ 1 . (23.3) 



The left hand side of (23.3) is a meromorphic function, with the leading zero 
at z = |A|. So what? 



Example 23.3 Meromorphic functions and exponential convergence: As an 

illustration of how exponential convergence of a truncated series is related to analytic 
properties of functions, consider, as the simplest possible example of a meromorphic 
function, the ratio 

z — a 
Kz) = -— 
z-b 

with a, b real and positive and a < b. Within the spectral radius |z| < b the function h 
can be represented by the power series 



Kz) = 



o-kzl'.. 



k=0 



where o-q - a/b, and the higher order coefficients are given by a-j - (a - b)/b-'^^ 
Consider now the truncation of order N of the power series 



a-kZ = - + 



z{a-b){\ -z'^/b'') 



,^0 ^ bHl-z/b) 

Let zn be the solution of the truncated series liNizN) = 0. To estimate the distance 
between a andzN it is sufficient to calculate liN{a). It is of order {a/b)'^^\ so finite order 
estimates converge exponentially to the asymptotic value. 



This example shows that: (1) an estimate of the leading pole (the leading 
eigenvalue of £.) from a finite truncation of a trace formula converges exponen- 
tially, and (2) the non-leading eigenvalues of £, lie outside of the radius of con- 
vergence of the trace formula and cannot be computed by means of such cycle 
expansion. However, as we shall now see, the whole spectrum is reachable at no 
extra effort, by computing it from a determinant rather than a trace. 



converg - 9nov2008 



ChaosBook.org versionl4.4.I, Apr 21 2013 



CHAPTER 23. WHY DOES IT WORK? 



455 



Figure 23.1: The Bernoulli shift map. 




Example 23.4 The spectral determinant for a single fixed point: 

determinant (19.3) follows from the trace formulas of example 23.2: 

CO / \ CO 

*'"-=^' = n('-siv) = Z<-')"e- 

where the cummulants Q„ are given explicitly by the Euler formula 

Qn = 



The spectral 



(23.4) 



1 - A-i 1 - A- 



1 - A-" 



exercise 23.3 



(23.5) 



The main lesson to glean from this simple example is that the cummulants Qn 
decay asymptotically /fli'ter than exponentially, as \-'>("-^)/^_ For example, if we 
approximate series such as (23.4) by the first 10 terms, the en^or in the estimate of 
the leading zero is w 1 / ! 

So far all is well for a rather boring example, a dynamical system with a single 
repelling fixed point. What about chaos? Systems where the number of unstable 
cycles increases exponentially with their length? We now turn to the simplest 
example of a dynamical system with an infinity of unstable periodic orbits. 



Example 23.5 Eigenfunction of Bernoulli shift map. (continued from example 11.7) The 
Bernoulli shift map figure 23.1 



fo(x) = 2x , 



xe/o = [0,1/2) 



/i(x) = 2jt:-l, xe/i=(l/2, 1] 



(23.6) 



models the 50-50% probability of a coin toss. The associated Perron-Frobenius oper- 
ator (16.9) assembles p{y) from its two prei mages 



(23.7) 



For this simple example the eigenfunctions can be written down explicitly: they coincide, 
up to constant prefactors, with the Bernoulli polynomials B„{x). These polynomials are 
generated by the Taylor expansion of the generating function 

&ix,t) = ^— = y Bkix)- , Boix) = 1 , Biix) ^x--,... 
e' - 1 j—i k\ 2 

k=i) 
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The Perron-Frobenius operator (23.7) acts on the generating function @ as 

^ ' k=\ 

hence each Bi,(x) is an eigenfunction of £ with eigenvaiue 1/2''. 

The fuii operator has two components corresponding to the two branches. For 
the n times iterated operator we have a fuii binary shift, and for each of the 2" branches 
the above caicuiations carry over, yieiding the same trace (2" - 1)^' for every cycie on 
iength n. Without further ado we substitute everything back and obtain the determinant. 



det 



(l-zi:) = exp[-X|5;^) = n(l-^)' (23.8) 

V 11=1 / k=0 



verifying that the Bernouiii poiynomiais are eigenfunctions with eigenvaiues 1,1/2,..., 
1/2", .... 



The Bernoulli map spectrum looks reminiscent of the single fixed-point spec- 
trum (23.2), with the difference that the leading eigenvalue here is 1, rather than 
1/|A|. The difference is significant: the single fixed-point map is a repeller, with 
escape rate (1.7) given by the X leading eigenvalue y = In |A|, while there is no 
escape in the case of the Bernoulli map. As already noted in discussion of the 
relation (19.23), for bounded systems the local expansion rate (here In |A| = In 2) section 19.4 
is balanced by the entropy (here In 2, the log of the number of preimages F^), 
yielding zero escape rate. 

So far we have demonstrated that our periodic orbit formulas are correct for 
two piecewise Unear maps in 1 dimension, one with a single fixed point, and one 
with a full binary shift chaotic dynamics. For a single fixed point, eigenfunctions 
are monomials in x. For the chaotic example, they are orthogonal polynomials on 
the unit interval. What about higher dimensions? We check our formulas on a 
2-dimensional hyperbolic map next. 



Example 23.6 The simplest of 2-dimensional maps - a single hyperbolic fixed 
point: We start by considering a very simpie iinear hyperboiic map with a singie 
hyperboiic fixed point, 

f(x)^(fi(xi,X2),f2(xi,X2))^(AsXi,A„X2), 0<|A,|<1, |A„| > 1 . 

The Perron-Frobenius operator (16.10) acts on the 2-dimensional density functions as 

Xpixuxi) = --^^p(xi/As,X2/A„) (23.9) 

What are good eigenfunctions? Cribbing the l-dimensionai eigenfunctions for the sta- 
ble, contracting xi direction from example 23. 1 is not a good idea, as under the iter- 
ation of X. the high terms in a Taylor expansion of p{xuX2) in the xi variable would 
get multiplied by exponentially exploding eigenvalues 1/A*. This makes sense, as in 
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the contracting directions liyperbolic dynamics crunclies up initiai densities, instead of 
smootliing tliem. So we guess instead tliat tlie eigenfunctions are of form 

(fk,k2(xi,X2) ^ xl^/x\'^\ ^1,^2=0,1,2,..., (23.10) 

a mixture of the Laurent series in the contraction x\ direction, and the Tayior series in 
the expanding direction, the X2 variabie. The action of Perron-Frobenius operator on 
this set of basis functions 

cr A^'' 

is smoothing, with the higher ki, k2 eigenvectors decaying exponentiaiiy faster, by 
^i'/Af,-^' factor in the eigenvaiue. One verifies by an expiicit caicuiation (undoing 
the geometric series expansions to iead to (19.9)) that the trace of -C indeed equais 
l/|det(l -M)| - 1/1(1 - A„)(l - A.s)| , from which it foiiows thataii our trace andspectrai 
determinant formuias appiy. The argument appiies to any hyperboiic map iinearized 
around the fixed point of form f{x\...., Xd) - (Aixi , A2X2, . . . , AdXd). 



So far we have checked the trace and spectral determinant formulas derived 
heuristically in chapters 18 and 19, but only for the case of 1 -dimensional and 
2-dimensional linear maps. But for infinite-dimensional vector spaces this game 
is fraught with dangers, and we have already been mislead by piecewise linear 
examples into spectral confusions: contrast the spectra of example 16.1 and ex- 
ample 17.4 with the spectrum computed in example 18.2. 

We show next that the above results do cany over to a sizable class of piece- 
wise analytic expanding maps. 



23.2 Evolution operator in a matrix representation 

The standard, and for numerical purposes sometimes very effective way to look at 
operators is through their matrix representations. Evolution operators are moving 
density functions defined over some state space, and as in general we can imple- 
ment this only numerically, the temptation is to discretize the state space as in 
sect. 16.3. The problem with such state space discretization approaches that they 
sometimes yield plainly wrong spectra (compare example 17.4 with the result of 
example 18.2), so we have to think through cai^efuUy what is it that we really 
measure. 

An expanding map f{x) takes an initial smooth density ^„{x), defined on a 
subinterval, stretches it out and overlays it over a larger interval, resulting in a new, 
smoother density 0„+i(x). Repetition of this process smoothes the initial density, 
so it is natural to represent densities (pn{x) by their Taylor series. Expanding 



k=Q f=0 
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and substitute the two Taylor series into (16.6): 



(j) = Cv) = I <ix S{y - fix)) (pnix) . 
JM 



The matrix elements follow by evaluating the integral 



v=o 



(23.11) 



we obtain a matrix representation of the evolution operator 



dx£(y,x)—=2_^—Ltk, k,k' =0,1,2,... 

k' 



which maps the x'^ component of the density of trajectories (pnix) into the y^' com- 
ponent of the density (pn+xiy) one time step later, with y = f{x). 

We already have some practice with evaluating derivatives 6^^\y) = -^^iy) from 
sect. 16.2. This yields a representation of the evolution operator centered on the 
fixed point, evaluated recursively in terms of derivatives of the map /: 



dx5^^\x- fix)) 



(23.12) 



x=f(x) 



1 



1 



\f'\\dxf'{x)j k\ 



'x^ 



The matrix elements vanish for { < k, so h is a lower triangular matrix. The 
diagonal and the successive off-diagonal matrix elements are easily evaluated it- 
eratively by computer algebra 



1 



lAIA'-- ' 



^k+l,k 



ik + l)\f" 
"2/t!|A|A^+2 ' 



For chaotic systems the map is expanding, |A| > 1. Hence the diagonal terms drop 
off exponentially, as 1/|A|^+^ , the terms below the diagonal fall off even faster, and 
truncating L to a finite matrix introduces only exponentially small enws. 

The trace formula (23.3) takes now a matrix form 



tr ■ 



l-z£ 



tr- 



1 -zL 



(23.13) 



converg - 9nov2008 



ChaosBook.org versionl4.4.I, Apr 21 2013 



CHAPTER 23. WHY DOES IT WORK? 



459 



Figure 23.2: A nonlinear one-branch repeller with a 
single fixed point w,. 



f(w) 
0.5 




In order to illustrate how this works, we work out a few examples. 

In example 23.7 we show that these results carry over to any analytic single- 
branch 1 -dimensional repeller. Further examples motivate the steps that lead to 
a proof that spectral determinants for general analytic 1 -dimensional expanding 
maps, and - in sect. 23.5, for 1-dimensional hyperbolic mappings - are also entire 
functions. 



Example 23.7 Perron-Frobenius operator in a matrix representation: As in ex- 
ample 23.1, we start with a map witli a single fixed point, but this time with a nonlin- 
ear piecewise analytic map f with a nonlinear inverse F = / sign of the derivative 
cr - cr{F') = F' /\F'\ , and the Perron-Frobenius operator acting on densities analytic in 
an open domain enclosing the fixed point x-Wq, 

Jl^iy) = Jdxdiy- fix)) cf>{x) = cr F'(y) cl>{F(y)) . 

Assume that F is a contraction of the unit disk in the complex plane, i.e., 

\F(z)\ <e<l and \F'{z)\ < C < oo for \z\ < I , 
and expand ^ in a polynomial basis with the Cauchy integral formula 

dw (p{w) r dw (f>(w) 



(23.14) 



0(z) = Xz'>„ = d)5- 



;i=0 



2ni w - z 



2ni w" 



Combining this with (23.22), we see that in this basis Perron-Frobenius operator £ is 
represented by the matrix 



r^r ^ V mr ^ r X (t F'{w){F(w))" 



(23.15) 



Taking the trace and summing we get: 

dw cr F'{w) 

n>{) 



2ni w - F{w) 



This integral has but one simple pole at the unique fixed point w* = F{w*) = /(w*). 
Hence exercise 23.6 



tr £ 



1 



cr F'jw*) _ 

1 - F'(w*) " |/'(w*)-l| 
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This super-exponential decay of cummulants ensures that for a repeller 
consisting of a single repelling point the spectral determinant (23.4) is entire in 
the complex z plane. 

In retrospect, the matrix representation method for solving the density evolu- 
tion problems is eminently sensible — after all, that is the way one solves a close 
relative to classical density evolution equations, the Schrodinger equation. When 
available, matrix representations for X enable us to compute many more orders 
of cumulant expansions of spectral determinants and many more eigenvalues of 
evolution operators than the cycle expensions approach. 

Now, if the spectral determinant is entire, formulas such as (19.25) imply that 
the dynamical zeta function is a meromorphic function. The practical import of 
this observation is that it guarantees that finite order estimates of zeroes of dyn- 
amical zeta functions and spectral determinants converge exponentially, or - in 
cases such as (23.4) - super-exponentially to the exact values, and so the cycle 
expansions to be discussed in chapter 20 represent a true perturbative approach to 
chaotic dynamics. 

Before turning to specifics we summarize a few facts about classical theory 
of integral equations, something you might prefer to skip on first reading. The 
purpose of this exercise is to understand that the Fredholm theory, a theory that 
works so well for the Hilbert spaces of quantum mechanics does not necessarily 
work for deterministic dynamics - the ergodic theory is much harder. 

fast track: 
sect. 23.4, p. 462 

23.3 Classical Fredholm theory 

He who would valiant be 'gainst all disaster 
Let him in constancy follow the Master. 
— John Bunyan, Pilgrim 's Progress 

The Perron-Frobenius operator 

£(Pix) = J dyS{x-f(y))m 

has the same appearance as a classical Fredholm integral operator 

'Kip(x) = f dy'K{x,y)<p(y), (23.16) 
Jm 

and one is tempted to resort to classical Fredholm theory in order to establish 
analyticity properties of spectral determinants. This path to enlightenment is 
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blocked by the singular nature of the kernel, which is a distribution, whereas the 
standard theory of integral equations usually concerns itself with regular kernels 
'7C(;c,j) e L^{Af). Here we briefly recall some steps of Fredholm theory, before 
working out the example of example 23.5. 

The general form of Fredholm integral equations of the second kind is 

ip(x) = f dynx,y)<p(y) + ^(x) (23.17) 
Jm 

where ^{x) is a given function in L^(M) and the kernel ^{x, y) e L^(At^) (Hilbert- 
Schmidt condition). The natural object to study is then the linear integral operator 
(23. 16), acting on the Hilbert space L^(A1): the fundamental property that follows 
from the L?iQ) nature of the kernel is that such an operator is compact, that is 
close to a finite rank operator.A compact operator has the property that for every 
6 > only a. finite number of linearly independent eigenvectors exist correspond- 
ing to eigenvalues whose absolute value exceeds 6, so we immediately realize 
(figure 23.5) that much work is needed to bring Perron-Frobenius operators into 
this picture. 

We rewrite (23.17) in the form 

Tip = ^, T = l-K. (23.18) 

The Fredholm alternative is now applied to this situation as follows: the equation 
Tip - ^ has a unique solution for every ^ £ L^{M) or there exists a non-zero 
solution of Tipo = 0, with an eigenvector of TC corresponding to the eigenvalue 1. 
The theory remains the same if instead of T we consider the operator Ta = l-A'K 
with /i 0. As 7C is a compact operator there is at most a denumerable set of A for 
which the second part of the Fredholm alternative holds: apart from this set the 
inverse operator (l-ATy^ exists and is bounded (in the operator sense). When A 
is sufficiently small we may look for a perturbative expression for such an inverse, 
as a geometric series 

(l-A'Kr^ = l + A'K + A^'K^ + ---= l + A'W, (23.19) 

where TC" is a compact integral operator with kernel 

"K^ix, y) ^ f dzi ... dzn-i 'K{x, zi) • • • %{zn-x, y) , 

JM"-i 

and is also compact, as it is given by the convergent sum of compact operators. 
The problem with (23.19) is that the series has a finite radius of convergence, 
while apart from a denumerable set of /I's the inverse operator is well defined. 
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A fundamental result in the theory of integral equations consists in rewriting the 
resolving kernel "W as a ratio of two analytic functions of A 



!Dix,y;A) 
D{A) 



If we introduce the notation 



X\... Xn 
yi ■ ■ -Jn 



nxuyi) ... nxuyn) 
'Kix„,yi) ... 'Kix„,y„) 



we may write the explicit expressions 



D{A) 



1 + 



exp 



2)ix,y;A) = 



TCI ^ 1+ 



1" r 

-V -trTC' 

m 



dzi ...dZn'K\ 



V m= 1 



(23.20) 



dzi...dZn'K\ I I' 



Zn 
Zn 



The quantity D{A) is known as the Fredholm determinant (see (19.24)):it is an 
entire analytic function of A, and D{A) = if and only if 1 //I is an eigenvalue of 
'K. 

Worth emphasizing again: the Fredholm theory is based on the compactness 
of the integral operator, i.e., on the functional properties (summability) of its ker- 
nel. As the Pen^on-Frobenius operator is not compact, there is a bit of wishful 
thinking involved here. 



23.4 Analyticity of spectral determinants 



They savored the strange warm glow of being much more 
ignorant than ordinary people, who were only ignorant of 
ordinary things. 

— Terry Pratchett 

Spaces of functions integrable L^, or square-integrable on interval [0, 1] 
are mapped into themselves by the Perron-Frobenius operator, and in both cases 
the constant function 0o = 1 is an eigenf unction with eigenvalue 1. If we focus 
our attention on we also have a family of eigenf unctions, 

4>e(y) = y,exp(27r/^3')-^ (23.21) 
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with complex eigenvalue 2 , parameterized by complex G with Re 9 > 0. By 
varying 9 one realizes that such eigenvalues fill out the entire unit disk. Such 
essential spectrum, the case ^ = of figure 23.5, hides all fine details of the 
spectrum. 

What's going on? Spaces and contain arbitrarily ugly functions, allow- 
ing any singularity as long as it is (square) integrable - and there is no way that 
expanding dynamics can smooth a kinky function with a non-differentiable singu- 
larity, let's say a discontinuous step, and that is why the eigenspectrum is dense 
rather than discrete. Mathematicians love to wallow in this kind of muck, but there 
is no way to prepare a nowhere differentiable initial density in a laboratory. The 
only thing we can prepare and measure are piecewise smooth (real-analytic) den- 
sity functions. 

For a bounded linear operator J?! on a Banach space il, the spectral radius 
is the smallest positive number pspec such that the spectrum is inside the disk of 
radius Pspec, while the essential spectral radius is the smallest positive number 
Pess such that outside the disk of radius p^ss the spectrum consists only of isolated 
eigenvalues of finite multiplicity (see figure 23.5). exercise 23.5 

We may shrink the essential spectrum by letting the Perron-Frobenius oper- 
ator act on a space of smoother functions, exactly as in the one-branch repeller 
case of sect. 23.1. We thus consider a smaller space, C^^", the space of k times 
differentiable functions whose k'th derivatives ai^e Holder continuous with an ex- 
ponent < a < 1 : the expansion property guarantees that such a space is mapped 
into itself by the Perron-Frobenius operator. In the strip O<Re0<^-i-a most (pg 
will cease to be eigenfunctions in the space C*^^"; the function 0„ survives only for 
integer valued 9 = n. In this way we anive at a finite set of isolated eigenvalues 
1, 2~', • • • , 2"*^, and an essential spectral radius p^^ = 2~'^'^'*'"K 

We follow a simpler path and restrict the function space even further, namely 
to a space of analytic functions, i.e., functions for which the Taylor expansion is 
convergent at each point of the interval [0, 1]. With this choice things turn out easy 
and elegant. To be more specific, let be a holomorphic and bounded function on 
the disk D = B(0, R) of radius R> centered at the origin. Our Perron-Frobenius 
operator preserves the space of such functions provided (1 -H /?)/2 < /? so all we 
need is to choose R > I. If Fg, s e {0,1}, denotes the s inverse branch of the 
Bernoulli shift (23.6), the corresponding action of the PeiTon-Frobenius operator 
is given by Jlshiy) - cr F'^iy) h o Fs(y), using the Cauchy integral formula along 
the dD boundary contour: 




(23.22) 



For reasons that will be made clear later we have introduced a sign cr = ± 1 of the 
given real branch |f (y)| = cr F'(y). For both branches of the Bernoulli shift s = I, 
but in general one is not allowed to take absolute values as this could destroy 
analyticity. In the above formula one may also replace the domain D by any 
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domain containing [0, 1] such that the inverse branches maps the closure of D into 
the interior of D. Why? simply because the kernel remains non-singular under 
this condition, i.e., w - F(y) + whenever w € dD and y € CI D. The problem 
is now reduced to the standard theory for Fredholm determinants, sect. 23.3. The 
integral kernel is no longer singular, traces and determinants are well-defined, and 
we can evaluate the trace of Xf by means of the Cauchy contour integral formula: 

£ dw <tF'(w) 



Ini w — F(w) 

Elementary complex analysis shows that since F maps the closure of D into its 
own interior, F has a unique (real-valued) fixed point x* with a multiplier strictly 
smaller than one in absolute value. Residue calculus therefore yields exercise 23.6 

crF'ix") 1 



l-F'(x*) \f'{x*)-l\' 
justifying our previous ad hoc calculations of traces using Dirac delta functions. 



Example 23.8 Perron-Frobenius operator in a matrix representation: As in ex- 
ample 23.1, we start with a map with a single fixed point, but this time with a nonlin- 
ear piecewise analytic map f with a nonlinear inverse F = sign of the derivative 

o- = a-{F') = F'I\F'\ 

mz)^ J dx6(z-f(x))<p{x)^CTF'(z)<f>(F{z)). 

Assume that F is a contraction of the unit d/sk, i.e., 

\F(z)\ <e<l and \F'(z)\ < C < oo for |z| < 1 , (23.23) 
and expand ^ in a polynomial basis by means of the Cauchy formula 

Combining this with (23.22), we see that in this basis X. is represented by the matrix 



CM mj ^ T cr F'(w){F(w)r 



(23.24) 



Taking the trace and summing we get: 

„ £ dw (J F'(w) 

tr i:=yL„„ = (b— ^ 

^—i J 2m w - r (w) 

;i>0 ^ ^ ' 

This integral has but one simple pole at the unique fixed point w* = F{w*) = /(w*). 
Henr.e 



Hence 

(T F'(w') ^ 1 
\-F'{w*) |/'(w*)-l| 
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We worked out a very specific example, yet our conclusions can be gener- 
alized, provided a number of restrictive requirements ai^e met by the dynamical 
system under investigation: exercise 23.6 

1) the evolution operator is multiplicative along the flow, 

2) the symbohc dynamics is a. finite subshift, 

3) all cycle eigenvalues are hyperbolic (exponentially bounded in 
magnitude away from 1), 

4) the map (or the flow) is real analytic, i.e., it has a piecewise ana- 
lytic continuation to a complex extension of the state space. 

These assumptions are romantic expectations not satisfied by the dynamical 
systems that we actually desire to understand. Still, they are not devoid of physical 
interest; for example, nice repeUers like our 3-disk game of pinball do satisfy the 
above requirements. 

Properties 1 and 2 enable us to represent the evolution operator as a finite 
matrix in an appropriate basis; properties 3 and 4 enable us to bound the size 
of the matrix elements and control the eigenvalues. To see what can go wrong, 
consider the following examples: 

Property 1 is violated for flows in 3 or more dimensions by the following 
weighted evolution operator 

£\y,x) = \A\x)f6(y-f{x)) , 

where A}{x) is an eigenvalue of the Jacobian matrix transverse to the flow. Semi- 
classical quantum mechanics suggest operators of this form with /? = l/2.The 
problem with such operators arises from the fact that when considering the Ja- 
cobian matrices Jab = JcJh for two successive trajectory segments a and b, the 
corresponding eigenvalues are in general not multiplicative, Kab + ^a^b (unless 
a, b are iterates of the same prime cycle p, so JaJb = Jp'^''')- Consequently, this 
evolution operator is not multiplicative along the trajectory. The theorems require 
that the evolution be represented as a matrix in an appropriate polynomial basis, 
and thus cannot be applied to non-multiplicative kernels, i.e., kernels that do not 
satisfy the semi-group property £! = ^* . 

Property 2 is violated by the 1-dimensional tent map (see figure 23.3 (a)) 

f{x) - a{\ - |1 - 2x\) , 1/2 < a < 1 . 

All cycle eigenvalues are hyperbolic, but in general the critical point Xc = 1/2 
is not a pre-periodic point, so there is no finite Markov partition and the sym- 
bolic dynamics does not have a finite grammar (see sect. 12.4 for definitions). In 
practice, this means that while the leading eigenvalue of £, might be computable. 
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Figure 23.3: (a) A (hyperbolic) tent map without 
a finite Maitov partition, (b) A Markov map with 
a marginal fixed point. 




the rest of the spectrum is very hard to control; as the parameter a is varied, the 
non-leading zeros of the spectral determinant move wildly about. 

Property 3 is violated by the map (see figure 23.3 (b)) 




X + 2x^ , X e Iq = [0, 5] 
2-2x , xeh ^[^,1] 



Here the interval [0, 1] has a Markov partition into two subintervals Iq and /i, and 
/ is monotone on each. However, the fixed point at x = has marginal stability 
Ao = 1, and violates condition 3. This type of map is called "intermittent" and 
necessitates much extra work. The problem is that the dynamics in the neighbor- 
hood of a marginal fixed point is very slow, with correlations decaying as power 
laws rather than exponentially. We will discuss such flows in chapter 24. 

Property 4 is required as the heuristic approach of chapter 18 faces two major 
hurdles: 

1. The trace (18.8) is not well defined because the integral kernel is singular. 

2. The existence and properties of eigenvalues are by no means clear. 

Actually, property 4 is quite restrictive, but we need it in the present approach, 
so that the Banach space of analytic functions in a disk is preserved by the Perron- 
Frobenius operator. 

In attempting to generalize the results, we encounter several problems. First, 
in higher dimensions life is not as simple. Multi-dimensional residue calculus is 
at our disposal but in general requires that we find poly-domains (direct product 
of domains in each coordinate) and this need not be the case. Second, and per- 
haps somewhat surprisingly, the 'counting of periodic orbits' presents a difficult 
problem. For example, instead of the Bernoulli shift consider the doubling map 
(11.8) of the circle, x 2x mod I, x e R/Z. Compared to the shift on the interval 
[0, 1] the only difference is that the endpoints and 1 are now glued together. Be- 
cause these endpoints are fixed points of the map, the number of cycles of length 
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n decreases by 1. The determinant becomes: 



( X— 1 7" 2" — ]\ 
de,(l-zX) = exp-;^- — 

V n=l / 



= 1 - z. 



(23.25) 



The value z = 1 still comes from the constant eigenfunction, but the Bernoulli 
polynomials no longer contribute to the spectrum (as they are not periodic). Proofs 
of these facts, however, are difficult if one sticks to the space of analytic functions. 

Third, our Cauchy formulas a priori work only when considering purely ex- 
panding maps. When stable and unstable directions co-exist we have to resort to 
stranger function spaces, as shown in the next section. 



Proceeding to hyperbolic systems, one faces the following paradox: If / is an 
area-preserving hyperbolic and real-analytic map of, for example, a 2-dimensional 
torus then the Perron-Frobenius operator is unitary on the space of functions, 
and its spectrum is confined to the unit circle. On the other hand, when we 
compute determinants we find eigenvalues scattered around inside the unit disk. 
Thinking back to the Bernoulli shift example 23.5 one would like to imagine 
these eigenvalues as popping up from the spectrum by shrinking the function 
space. Shrinking the space, however, can only make the spectrum smaller so this 
is obviously not what happens. Instead one needs to introduce a 'mixed' function 
space where in the unstable direction one resorts to analytic functions, as before, 
but in the stable direction one instead considers a 'dual space' of distributions on 
analytic functions. Such a space is neither included in nor includes and we 
have thus resolved the paradox. However, it still remains to be seen how traces 
and determinants are calculated. 

The linear hyperbolic fixed point example 23.6 is somewhat misleading, as we 
have made explicit use of a map that acts independently along the stable and unsta- 
ble directions. For a more general hyperbolic map, there is no way to implement 
such direct product structure, and the whole argument falls apart. Her comes an 
idea; use the analyticity of the map to rewrite the Perron-Frobenius operator acting 
as follows (where cr denotes the sign of the derivative in the unstable direction): 



23.5 Hyperbolic maps 



I can give you a definion of a Banach space, but I do not 
know what that means. 

— Federico Bonnetto, Banach space 



(H.H. Rugh) 



-Ch{Zi,Z2) = 




cr /i(wi, W2) dw[ dw2 



(23.26) 



(zi -/i(wi,W2)(/2(wi,W2) -Z2) 2ni 2ni 
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Figure 23.4: For an analytic hyperbolic map, specify- 
ing the contracting coordinate wi, at the initial rectangle 
and the expanding coordinate z,, at the image rectangle 
defines a unique trajectory between the two rectangles. 
In particular, Wy and zi, (not shown) are uniquely spec- 
ified. 



Here the function (p should belong to a space of functions analytic respectively 
outside a disk and inside a disk in the first and the second coordinates; with the 
additional property that the function decays to zero as the first coordinate tends 
to infinity. The contour integrals are along the boundaries of these disks. It is 
an exercise in multi-dimensional residue calculus to verify that for the above lin- 
ear example this expression reduces to (23.9). Such operators form the building 
blocks in the calculation of traces and determinants. One can prove the following: 

Theorem: The spectral determinant for 2-dimensional hyperbolic analytic maps 

is entire. remark 23.8 

The proof, apart from the Markov property that is the same as for the purely 
expanding case, relies heavily on the analyticity of the map in the explicit con- 
struction of the function space. The idea is to view the hyperbolicity as a cross 
product of a contracting map in forward time and another contracting map in back- 
ward time. In this case the Markov property introduced above has to be elaborated 
a bit. Instead of dividing the state space into intervals, one divides it into rectan- 
gles. The rectangles should be viewed as a direct product of intervals (say hori- 
zontal and vertical), such that the forward map is contracting in, for example, the 
horizontal direction, while the inverse map is contracting in the vertical direction. 
For Axiom A systems (see remark 23.8) one may choose coordinate axes close 
to the stable/unstable manifolds of the map. With the state space divided into 
A'^ rectangles {Ali, ■ ■ ■ ,Mn}, Mi = x 71' one needs a complex extension 
D'' X DV, with which the hyperbolicity condition (which simultaneously guaran- 
tees the Markov property) can be formulated as follows: 

Analytic hyperbolic property: Either /(A1,) n Int(Aly) = 0, or for each pair 
Wh £ C1(D|'), Zv £ Cl(Dp there exist unique analytic functions of Wh,Zv' Wv = 
Wviwh,Zv) e Int(Dp, Zh = Zh(wh,Zv) e Int(D^O, such that /(w/,, w,,) = (Zh,Zv)- 
Furthermore, if w/j £ if and Zv £ Tj, then w,, £ IJ and Zh £ I'j (see figure 23.4). 

In plain English, this means for the iterated map that one replaces the coor- 
dinates Zh,Zv at time n by the contracting pair Zh, Wy, where Wy is the contracting 
coordinate at time n + I for the 'partial' inverse map. 

In two dimensions the operator in (23.26) acts on functions analytic outside 
in the horizontal direction (and tending to zero at infinity) and inside DJ' in 

the vertical direction. The contour integrals are precisely along the boundaries of 

these domains. 
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A map / satisfying the above condition is called analytic hyperbolic and the 
theorem states that the associated spectral determinant is entire, and that the trace 
formula (18.8) is correct. 

Examples of analytic hyperbolic maps are provided by small analytic pertur- 
bations of the cat map, the 3-disk repeller, and the 2-dimensional baker's map. 



By now we appreciate that any honest attempt to look at the spectral prop- 
erties of the Perron-Frobenius operator involves hard mathematics, but the reward 
is of this effort is that we are able to control the analyticity properties of dynamical 
zeta functions and spectral determinants, and thus substantiate the claim that these 
objects provide a powerful and well-founded theory. 

Often (see chapter 17) physically important part of the spectrum is just the 
leading eigenvalue, which gives us the escape rate from a repeller, or, for a gen- 
eral evolution operator, formulas for expectation values of observables and their 
higher moments. Also the eigenfunction associated to the leading eigenvalue has 
a physical interpretation (see chapter 16): it is the density of the natural measures, 
with singular measures ruled out by the proper choice of the function space. This 
conclusion is in accord with the generalized PeiTon-Frobenius theorem for evolu- 
tion operators. In a finite dimensional setting, the statement is: remark 23.7 

• Perron-Frobenius theorem: Let be a non-negative matrix, such that 
some finite n exists for which any initial state has reached any other state, 
{L")ij > ViJ: then 

1. The maximal modulus eigenvalue is non-degenerate, real, and posi- 



2. The corresponding eigenvector (defined up to a constant) has non- 
negative coordinates. 

We may ask what physical information is contained in eigenvalues beyond the 
leading one: suppose that we have a probability conserving system (so that the 
dominant eigenvalue is 1), for which the essential spectral radius satisfies < 
Pess < ^ < 1 on some Banach space S. Denote by P the projection con^esponding 
to the part of the spectrum inside a disk of radius 6. We denote by Ai, A2 ... , Am 
the eigenvalues outside of this disk, ordered by the size of their absolute value, 
with A\ = I. Then we have the following decomposition 



23.6 Physics of eigenvalues and eigenfunctions 



five. 



M 




(23.27) 
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when Lj are (finite) matrices in Jordan canomical form (Lq = is a [1 x 1] matrix, 
as Aq is simple, due to the PeiTon-Frobenius theorem), whereas ipi is a row vector 
whose elements form a basis on the eigenspace corresponding to Ai, and ip* is 
a column vector of elements of S* (the dual space of linear functionals over S) 
spanning the eigenspace of £,* corresponding to Aj. For iterates of the Perron- 
Frobenius operator, (23.27) becomes 

M 

rv = ^"-A/^- + prv ■ (23.28) 

i=l 

If we now consider, for example, correlation between initial (p evolved n steps and 
final ^, 

{^\r\ip} = [ dy ^(y) (£» (y) = f dw(^o f")(wMw) , (23.29) 
JM JM 

it follows that 

L 

{^\£"\,p} = A"^cji(^, 'P) + Y, f) + 0(9") , (23.30) 

i=2 

where 

ojf\^,<p) = f dy^(y)ilfiLiri^. 
JM 

The eigenvalues beyond the leading one provide two pieces of information: 
they rule the convergence of expressions containing high powers of the evolution 
operator to leading order (the A\ contribution). Moreover if (jD\{^,(p) = then exercise 23.7 
(23.29) defines a coiTclation function: as each term in (23.30) vanishes exponen- 
tially in the n oo limit, the eigenvalues Ai, . . . , Am determine the exponential 
decay of correlations for our dynamical system. The prefactors co depend on the 
choice of functions, whereas the exponential decay rates (given by logarithms of 
Ai) do not: the correlation spectrum is thus a universal property of the dynamics 
(once we fix the overall functional space on which the Perron-Frobenius operator 
acts). 



Example 23.9 Bernoulli shift eigenfunctions: Let us revisit the Bernoulli shift ex- 
ample (23.6) on the space of analytic functions on a disk: apart from the origin we have 
only simple eigenvalues - 2^*, ^ = 0, 1, . . .. The eigenvalue Ao = 1 corresponds to 
probability conservation: the corresponding eigenfunction Bo{x) - 1 indicates that the 
natural measure has a constant density over the unit interval. If we now take any ana- 
lytic function ri(x) with zero average (with respect to the Lebesgue measure), it follows 
that 0)1(77,77) = 0, and from (23.30) the asymptotic decay of the correlation function is 
(unless also 011(77, 77) = 0) 

C,,,(7i) ~ exp(-72log2). (23.31) 
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Thus, -log/ti gives the exponential decay rate of correlations (with a prefactor that 
depends on the choice of the function). Actually the Bernoulli shift case may be treated 
exactly, as for analytic functions we can employ the Euler-t^acLaurin summation for- 
mula 

7j(z) = dwiiiw) + y ^ —B„Az). (23.32) 

Jo rn\ 

As we are considering functions with zero average, we have from (23.29) and the fact 
that Bernoulli polynomials are eigenvectors of the Perron-Frobenius operator that 

(2-"')"(/7'™'(l)-/7*"''(0)) ri 



C,,,,,{n) = 2j i f dzT]{z)B,„{z) . 



«7— 1 



The decomposition (23.32) is also useful in realizing that the linear functionals tff* are 
singular objects: if we write it as 



m=Q 

we see that these functionals are of the form 



Jo 

where 

>F,.(w) = ^-77— - 1) - (5<''"'*(w)) , (23.33) 

when i > 1 and^oiw) - 1. This representation is only meaningful when the function s 
is analytic in neighborhoods ofw,w-l. 



23.7 Troubles ahead 



The above discussion confirms that for a series of examples of increasing gener- 
ality formal manipulations with traces and determinants are justified: the Perron- 
Frobenius operator has isolated eigenvalues, the trace formulas are explicitly ver- 
ified, and the spectral determinant is an entire function whose zeroes yield the 
eigenvalues. Real life is harder, as we may appreciate through the following 
considerations: 



• Our discussion tacitly assumed something that is physically entirely reason- 
able: our evolution operator is acting on the space of analytic functions, i.e., 
we are allowed to represent the initial density p(x) by its Taylor expansions 
in the neighborhoods of periodic points. This is however far from being the exercise 23.1 
only possible choice: mathematicians often work with the function space 
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Figure 23.5: Spectrum of the Perron-Frobenius oper- 
ator acting on the space of C*"^" Holder-continuous 
functions: only k isolated eigenvalues remain between 
the spectral radius, and the essential spectral radius 
which bounds the "essential," continuous spectrum. 




spWtral radius j.^iated eige/value 



£k+a^ i.e., the space of k times differentiable functions whose ^'th deriva- 
tives are Holder continuous with an exponent < a < 1 : then every y'' with 
Re 7/ > A: is an eigenfunction of the Perron-Frobenius operator and we have 

This spectrum differs markedly from the analytic case: only a small number 
of isolated eigenvalues remain, enclosed between the spectral radius and a 
smaller disk of radius 1/|A|^''^\ see figure 23.5. In literature the radius of 
this disk is called the essential spectral radius. 

In sect. 23.4 we discussed this point further, with the aid of a less trivial 
1 -dimensional example. The physical point of view is complementary to 
the standard setting of ergodic theory, where many chaotic properties of a 
dynamical system are encoded by the presence of a continuous spectrum, 
used to prove asymptotic decay of correlations in the space of L} square- 
integrable functions. exercise 23.2 

A deceptively innocent assumption is hidden beneath much that was dis- 
cussed so far: that (23.1) maps a given function space into itself. The ex- 
panding property of the map guarantees that: if f{x) is smooth in a do- 
main D then f{x/A) is smooth on a larger domain, provided | A| > 1 . For 
higher-dimensional hyperbolic flows this is not the case, and, as we saw in 
sect. 23.5, extensions of the results obtained for expanding 1-dimensional 
maps are highly nontrivial. 

It is not at all clear that the above analysis of a simple one-branch, one fixed 
point repeller can be extended to dynamical systems with Cantor sets of 
periodic points: we showed this in sect. 23.4. 



Resume 



Examples of analytic eigenf unctions for 1-dimensional maps are seductive, and 
make the problem of evaluating ergodic averages appear easy; just integrate over 
the desired observable weighted by the natural measure, right? No, generic natural 
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measure sits on a fractal set and is singular everywhere. The point of this book 
is that you never need to construct the natural measure, cycle expansions will do 
that job. 

A theory of evaluation of dynamical averages by means of trace formulas 
and spectral determinants requires a deep understanding of their analyticity and 
convergence. We worked here through a series of examples: 

1 . exact spectrum (but for a single fixed point of a linear- map) 

2. exact spectrum for a locally analytic map, matrix representation 

3. rigorous proof of existence of discrete spectrum for 2-dimensional hyper- 
bolic maps 

In the case of especially well-behaved "Axiom A" systems, where both the 
symbolic dynamics and hyperbolicity ai^e under control, it is possible to treat 
traces and determinants in a rigorous fashion, and strong results about the ana- 
lyticity properties of dynamical zeta functions and spectral determinants outlined 
above follow. 

Most systems of interest are not of the "axiom A" category; they are neither 
purely hyperbolic nor (as we have seen in chapters 11 and 12 ) do they have 
finite grammar. The importance of symbolic dynamics is generally grossly under 
appreciated; the crucial ingredient for nice analyticity properties of zeta functions 
is the existence of a finite grammar (coupled with uniform hyperbolicity). 

The dynamical systems which are really interesting - for example, smooth 
bounded Hamiltonian potentials - ai^e presumably never fully chaotic, and the 
central question remains: How do we attack this problem in a systematic and 
controllable fashion? 
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Theorem: Conjecture 3 with technical hypothesis is true 
in a lot of cases. 
— M. Shub 



Commentary 

Remark 23.1 Surveys of rigorous theory. We recommend the references hsted in re- 
mark 1 . 1 for an introduction to the mathematical literature on this subject. For a physicist, 
Driebe's monograph [1.19] might be the most accessible introduction into mathematics 
discussed briefly in this chapter There are a number of reviews of the mathematical ap- 
proach to dynamical zeta functions and spectral determinants, with pointers to the original 
references, such as refs. [23.1, 23.2]. An alternative approach to spectral properties of the 
Perron-Frobenius operator is given in ref. [23.3]. 

Ergodic theory, as presented by Sinai [23.14] and others, tempts one to describe the 
densities on which the evolution operator acts in terms of either integrable or square- 
integrable functions. For our purposes, as we have already seen, this space is not suitable. 
An introduction to ergodic theory is given by Sinai, Kornfeld and Fomin [23.15]; more 
advanced old-fashioned presentations are Walters [23.12] and Denker, Grillenberger and 
Sigmund [23.16]; and a more formal one is given by Peterson [23.17]. 

Remark 23.2 Fredholm theory. Our brief summary of Fredholm theory is based 
on the exposition of ref. [23.4]. A technical introduction of the theory from an operator 
point of view is given in ref. [23.5]. The theory is presented in a more general form in 
ref. [23.6]. 

Remark 23.3 Bernoulli shift. For a more in-depth discussion, consult chapter 3 
of ref. [1.19]. The extension of Fredholm theory to the case or Bernoulli shift on C^"^" 
(in which the Perron-Frobenius operator is not compact - technically it is only quasi- 
compact. That is, the essential spectral radius is strictly smaller than the spectral radius) 
has been given by Ruelle [23.7]: a concise and readable statement of the results is con- 
tained in ref. [23.8]. We see from (23.31) that for the Bernoulli shift the exponential 
decay rate of correlations coincides with the Lyapunov exponent: while such an identity 
holds for a number of systems, it is by no means a general result, and there exist explicit 
counterexample s . 

Remark 23.4 Hyperbolic dynamics. When dealing with hyperbolic systems one 
might try to reduce to the expanding case by projecting the dynamics along the unstable 
directions. As mentioned in the text this can be quite involved technically, as such unstable 
foliations are not characterized by strong smoothness properties. For such an approach, 
see ref. [23.3]. 

Remark 23.5 Spectral determinants for smooth flows. The flieorem on page 467 also 
applies to hyperbolic analytic maps in d dimensions and smooth hyperbolic analytic flows 
in{d + 1) dimensions, provided that the flow can be reduced to a piecewise analytic map 
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by a suspension on a Poincare section, complemented by an analytic "ceiling" function 
(3.5) that accounts for a variation in the section return times. For example, if we take 
as the ceiling function g{x) = e'^'**, where T{x) is the next Poincare section time for a 
trajectory staring at x, we reproduce the flow spectral determinant (19.13). Proofs are 
beyond the scope of this chapter 

Remark 23.6 Explicit diagonalization. For 1 -dimensional repellers a diagonalization 
of an explicit truncated L„,„ matrix evaluated in a judiciously chosen basis may yield many 
more eigenvalues than a cycle expansion (see refs. [23.10, 23.11]). The reasons why one 
persists in using periodic orbit theory are partially aesthetic and partially pragmatic. The 
explicit calculation of L,„„ demands an explicit choice of a basis and is thus non-invariant, 
in contrast to cycle expansions which utilize only the invariant information of the flow. In 
addition, we usually do not know how to construct L,„„ for a realistic high-dimensional 
flow, such as the hyperbolic 3-disk game of pinball flow of sect. 1.3, whereas periodic 
orbit theory is true in higher dimensions and straightforward to apply. 

Remark 23.7 Perron-Frobenius tlieorem. A proof of the Perron-Frobenius theo- 
rem may be found in ref. [23.12]. For positive transfer operators, this theorem has been 
generalized by Ruelle [23.13]. 

Remark 23.8 Axiom A systems. The proofs in sect. 23.5 follow the thesis 

work of H.H. Rugh [23.9, 23.18, 23.19]. For a mathematical introduction to the subject, 
consult the excellent review by V. Baladi [23.1]. It would take us too far afield to give 
and explain the definition of Axiom A systems (see refs. [1.27, 1.28]). Axiom A implies, 
however, the existence of a Markov partition of the state space from which the properties 
2 and 3 assumed on page 456 follow. 

Remark 23.9 Left eigenfunctions. We shafl never use an expUcit form of left eigen- 
functions, corresponding to highly singular kernels like (23.33). Many details have been 
elaborated in a number of papers, such as ref. [23.20], with a daring physical interpreta- 
tion. 

Remark 23.10 Ulam's idea. The approximation of Perron-Frobenius operator defined 
by (16.14) has been shown to reproduce the spectrum for expanding maps, once finer 
and finer Markov partitions are used [23.21]. The subtle point of choosing a state space 
partitioning for a "generic case" is discussed in ref. [23.22]. 
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Exercises 



23.1. What space does £ act on? Show that (23.2) is a 
complete basis on the space of analytic functions on a 
disk (and thus that we found the complete set of eigen- 
values). 

23.2. What space does X act on? What can be said about 
the spectrum of (23.1) on L'[0, 1]? Compare the result 
with figure 23.5. 

23.3. Euler formula. Derive the Euler formula (23.5), 
\u\ < 1: 

CO 9 

n,, t t U 

(1 + fM*) = 1 + + 

, „ ' l-u (1 -m)(1 -m2) 



ife=0 



(1 -m)(1 -m2)(1 -m3) 

CO t(l-l) 



(1 -m)---(1 -M*^) 



k=0 



23.4. 2-dimensional product expansion**. We conjecture 
that the expansion corresponding to exercise 23.3 is in 
the 2-dimensional case given by 



]^(l+f«^) 



S (1 -u 



Fkiu) 



k=0 

= 1 + 



(1 -m)^(1 -m2)2..-(1 -M*) 



,k\2 



I 



2m 



(l-u)- (I -«)2(l -m2)2 

u^(l +4u + u^) 



r 



(l -m)2(1 -m2)2(1 _„3)2 



k=0 



Fk{u) is a polynomial in m, and the coefficients fall off 
asymptotically as C„ ^ m" Verify; if you have a proof 
to all orders, e-mail it to the authors. (See also solu- 
tion 23.3). 

23.5. BernoulU shift on L spaces. Check that the family 
(23.21) belongs to L'([0, 1]). What can be said about 
the essential spectral radius on L^([0, 1])? A useful ref- 
erence is [23.24]. 

23.6. Cauchy integrals. Rework all complex analysis steps 
used in the BernoulU shift example on analytic functions 
on a disk. 

23.7. Escape rate. Consider the escape rate from a strange 
repeller: find a choice of trial functions ^ and ip such 
that (23.29) gives the fraction on particles surviving after 
n iterations, if their initial density distribution is pq{x). 
Discuss the behavior of such an expression in the long 
time limit. 
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